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EFFECTS  OF  CATHODE  SURFACE  ROUGHNESS 
ON  THE  QUALITY  OF  ELECTRON  BEAMS 


I.  Introduction 

A  critical  factor  in  the  operation  of  coherent  radiation  sources  such  as 
free  electron  lasers  (FEL),  gyrotrons,  etc.,  is  the  quality  of  the  electron 
beams. It  has  been  known  for  some  time  that  the  beam  quality  is  limited  by 
the  following  factors: the  electron  thermal  velocities  which  accompany  the 
finite  cathode  temperature,  the  roughness  on  the  cathode  surfaces,  patchiness 
and  non-uniform  emission,  asymmetries  and  non-uniformities  in  the  confining 
fields,  non-adiabatlc  fields,  space  charge  fields,  non-linear  forces,  etc. 
Among  these,  the  thermal  velocities  are  considered  unimportant  for  gyrotrons 
and  FELs  driven  by  induction  accelerators.  Non-uniform  and  non-adiabatic 
fields  may  be  considerably  reduced  with  careful  gun  design.  37  Space  charge 
effects  may  be  corrected  by  the  well-known  Pierce  method,  at  least  for  one  set 
of  parameters.  It  then  appears  that  the  problems  which  are  cathode-related 
such  as  patchiness,  non-uniform  emission,  roughness,  and  non-constant  work 
function,  etc.,  may  in  many  cases  place  a  fundamental  limit  on  the  brightness 


of  an  electron  beam. 


9.7 


In  this  paper,  we  seek  to  provide  some  preliminary 


assessment  of  the  effects  of  surface  roughness  on  the  beam  quality.  Analytic 
scaling  laws  are  derived  which  yield  the  upper  bound  on  the  beam 
emittance10-1^  in  terms  of  the  roughness  scales,  for  both  temperature  limited 
and  space  charge  limited  regimes.  These  scaling  laws  have  been  supported  by 
numerical  integration  of  the  electron  trajectories. 
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Velocity  spread  due  to  surface  roughness  has  been  studied  previously  for 
the  temperature  limited  regime.  Tsimring^  estimated  that  surface  roughness  of 
5  um  size  may  lead  to  a  5  per  cent  velocity  spread  in  magnetron  injection 
guns.  More  recently,  Kapetanakos^  estimated  that  a  fifty  micron  roughness  may 
result  in  a  normalized  emittance  of  0.25  cm-rad  for  a  cathode  of  5  cm  radius. 
These  estimates  pointed  to  the  significance  of  surface  roughness  on  the 
performance  of  an  electron  beam.  Both  Tsimring  and  Kapetanakos  conjectured 
that  the  deleterious  effects  of  surface  roughness  may  be  reduced  if  the 
cathode  is  operated  in  the  space  charge  limited  regime.  It  is  these  earlier 
speculations,  together  with  the  pressing  need  to  address  the  beam  quality 
issue,  which  prompted  the  more  detailed  study  to  be  reported  in  this  paper. 

Here,  we  shall  only  consider  non-immersed  cathodes,  where  the  external 
magnetic  field  is  absent.  Since  the  surface  roughness  has  a  typical  dimension 
on  the  order  of  tens  of  microns  (or  less),  it  is  anticipated  that  it3 
influence  on  the  electron  trajectories  will  be  limited  to  a  distance  on  the 
order  of  a  few  times  this  dimension,  over  which  the  self-magnetic  field  of  the 
electron  beam  may  certainly  be  ignored.  Thus,  magnetic  field  effects  are 
ignored  altogether  in  the  present  study  and  the  beam  tranverse  velocities  are 
caused  only  by  the  electrostatic  fields.  This  observation  provides 
considerable  simplification.  For  example,  it  allows  a  two-dimensional 
formulation.  As  a  result,  a  reliable,  yet  tractable  solution  in  the 
temperature  limited  regime  may  be  constructed. 

In  the  temperature  limited  regime,  the  space  charge  effects  are 
negligible.  The  electrostatic  field  is  simply  the  vacuum  field  which  may  be 
readily  solved  by  conformal  mapping  if  a  convenient,  two-dimensional  model  is 
used  to  represent  a  rough  surface.  The  electron  trajectories  subject  to  this 
vacuum  field  may  be  routinely  computed  numerically.  A  scaling  law  may  be 
regarded  as  reliable  if  it  agrees  with  numerical  computation  over  a  wide  range 


of  parameters. 


Matters  are  substantially  more  complicated  in  the  space  charge  limited 
regime,  where  the  space  charge  density  is  highest  near  the  cathode  surface. 

An  accurate  assessment  of  the  effects  of  surface  roughness  would  have 
required,  at  a  minimum,  a  two-dimensional  solution  of  the  Child-Langmuir  type 
over  a  bumpy  suface.  Such  a  solution  necessarily  reflects  the  self- 
consistency  between  charge  distribution  and  electron  trajectories,  and  an 
analytic  solution  does  not  seem  to  have  been  constructed.  It  is  not  the 
objective  of  this  paper  to  construct  such  a  solution.  Rather,  we  attempt  to 
provide  a  semi-quantitative  answer — based  on  what  we  believe  to  be  a 
reasonable  procedure — as  to  what  extent  the  space  charge  limited  operation 
would  alleviate  the  surface  roughness  problem.  Essentially,  we  used  the  one¬ 
dimensional  Child-Langmuir  field,  mapped  it  onto  the  bumpy  surface  as  the 
lowest  order  approximation,  and  then  calculate  the  electron  trajectories 
subject  to  this  transformed  field.  Note  that  this  transformed  field,  while  it 
is  not  consistent  with  the  trajectories,  does  preserve  certain  expected 
properties.  A  scaling  law  is  next  constructed  and  compared  with  numerical 
integration  based  on  the  procedure.  Agreement  is  observed.  But  this 
agreement  does  not  eliminate  the  possibility  that  our  analysis  of  the  space 
charge  limited  regime  may  be  of  limited  validity  so  that  the  resulting  scaling 
law  should  be  regarded  as  tentative.  However,  one  estimate  shows  that  the 
error  committed  is  on  the  order  of  thirty  per  cent. 

The  result  of  the  present  study  is  summarized  as  follows.  The  transverse 
momentum  caused  by  surface  roughness  is  computed  numerically  and  estimated 
analytically.  The  maximum  values  are  approximately  given  by  Eqs.  (15)  and 
(16),  respectively  for  the  temperature  limited  and  space  charge  limited 
regime.  These  formulas  have  been  favorably  tested  against  numerical 
computation  over  a  wide  range  of  parameters.  From  these  calculations,  we 
conclude  that  the  beam  emlttance,  due  to  surface  roughness  alone,  is  lower  for 


apace  charge  limited  operation  than  for  temperature  limited  operation  by  a 
factor  of  two  to  five. 

The  scaling  laws  (15),  (16)  provide  an  immediate  assessment  of  the 
effects  of  surface  roughness.  An  example  suffices  to  illustrate  the  order  of 
magnitude.  For  an  anode-cathode  voltage  drop  of  250  keV  over  a  distance  of 
5  cm,  say,  surface  roughness  of  5  uni  size  leads  to  an  equivalent  maximum 
transverse  temperature  of  10  eV  if  the  cathode  is  operated  in  the  temperature 
limited  regime,  but  of  0.55  eV  in  the  space  charge  limited  regime.  Since 
these  equivalent  temperatures  are  considerably  higher  than  the  normal 
operating  temperature  of  cathodes,  surface  roughness  (and  perhaps  cathode 
patchiness  and  non-uniform  emission)  deserves  attention  in  those  cases  where 
thermal  effects  are  considered  serious.  On  the  other  hand,  in  most  existing 
induction  accelerators,  surface  roughness  alone  is  unlikely  to  be  the  only 
deciding  factor  in  the  emittance  of  a  post-accelerated  electron  beam. 

In  Section  II,  we  describe  the  model  in  more  detail  and  present  the 
results  for  the  temperature  limited  regime.  In  Section  III,  the  space  charge 
limited  regime  is  examined.  The  analytic  scaling  laws  are  derived  in  Section 
IV,  and  the  concluding  remarks  will  be  given  in  the  last  section.  Some 
technical  details  are  outlined  in  the  Appendices. 


II.  The  Model  and  the  Temperature  Limited  Regime 

Consider  a  cathode  surface  whose  roughness  is  characterized  by  height  h, 
width  w  and  separation  L  [Fig.  1,  Top].  If  L  >>  h,  w,  the  bumps  on  the  rough 
surface  are  uncorrelated  and  we  can  consider  the  emission  of  a  single  bump 
[Fig  1,  Bottom].  For  simplicity,  we  consider  a  two-dimensional  model  of  the 
protusion,  whose  profile  is  described  by  the  equation 

y  -  b[l  +  - 2^/2  ‘  b*  (1  ) 

b  +  x 


Here,  a  and  b  are  constants  depending  on  the  height  h  and  w,  the  full  width  at 
half  maximum.  The  relation  between  a,  b  and  h,  w  is  given  in  Appendix  A.  The 
choice  of  Eq.  (1)  for  the  surface  is  motivated  by  the  simple  transformation 
f(z)  -  (z2  +  1)1/2  where  z  -  x  +  iy.  [See  Appendix  A], 

Of  interest  is  the  degree  of  divergence  of  an  electron  emitted  from  the 
rough  surface  described  in  Eq.  (1).  Assuming  that  the  anode  is  located  at 
y  ■  D,  where  D  »  h,  w,  the  vacuum  electric  field  Ex,  Ey  is  approximately 
given  by  the  expression  [See  Appendix  A] 

-iE  [x  -  i(y  +  b)] 

E  +  iE  -  - 5 - 2_ - - - -  (2) 

y  [(a2  -  (y  +  b)2  +  x2)  -  2ix(y  +  b)] 


Equation  (2)  is  valid  if  the  space  charge  effect  is  negligible.  Here,  E0  is 
the  average  electric  field  $/D,  $  is  the  anode-cathode  voltage  drop.  The 
instantaneous  position  (x,y)  and  the  normalized  instantaneous  momentum  (px, 
py)  are  governed  by  the  following  equations 


—  (Px) 

dt  lp./ 


_e_  ( 

m  c  v 


(3) 


dt  "y'  ~  '  / m  2  2.1/2  vp  ' 

y  (1  +  p  +  p  )  *y 

x  y 

If  an  electron  is  emitted  with  zero  velocity  the  initial  conditions  at  time 
t  -  0  are: 


x  -  xQ  <uQ) 

y  -  y0  (u0) 

Px  ■  o 
Py  -  0. 


In  these  equations,  e  is  the  electron  charge,  mQ  is  the  electron  rest  mass,  c 

is  the  speed  of  light,  c8  ,  c8  are  the  components  of  the  instantaneous 

x  y 

velocity  of  the  electron,  p  -  Y8  ,  p  -  Y8  with  Y  »  (1  +  p2  +  p2)1/2  being 

x  x  y  y  x  y 

the  relativistic  mass  factor,  and  uQ  in  Eq.  (5)  signifies  the  location  on  the 

surface  from  which  an  electron  is  emitted.  [See  Eq3.  (A3),  (A4)  of  Appendix  A 

Once  Eq,  h,  w  are  specified,  the  field  profile  is  determined  from  Eq.  (2) 

and  (x,y,px,py)  can  be  numerically  integrated  in  time  according  to  (3),  (4) 

until  the  electron  reaches  an  energy  equal  to  EqD  -  $.  Shown  in  Figs.  2,  3 

are  typical  numerical  results  for  various  electrons  as  they  reach  the  anode. 

The  transverse  displacement  (x  -  xQ)  as  a  function  of  uQ/uE  is  given  in  Fig.  2 

for  h  -  13  ym,  w  -  0.1  h,  E0  =  0.3  MeV/cra  and  various  values  of  <p.  Roughly, 

u  /u„  is  a  measure  of  the  distance  along  the  rough  surface  from  the  peak  of 
O  b 

the  bump.  More  precisely,  uQ  is  defined  in  Eq.  (A3,  A4)  in  terms  of  (xQ,  yQ) , 
the  initial  position  of  an  electron,  u0  =  0  corresponds  to  the  peak  of  the 
bump  and  uE  corresponds  to  a  point  at  the  tail  of  the  bump  at  which 
yo(uE)  »  n/10.  [See  appendix  A] 


The  data  may  be  understood  as  follows.  Higher  values  of  <j»  imply  a  larger 


value  of  Y  and  therefore  a  larger  value  of  YB  (x  -  x0)  as  shown  in  Fig.  2. 

Note,  however,  that  (x  -  xQ)  is  only  weakly  dependent  on  $.  Figure  3  shows 

the  value  Y88  evaluated  at  the  time  when  the  electrons  reach  the  anode. 

Here  0  ■  v  /v  .  At  u  -  0,  x*  »  0,  the  electric  field  has  no  transverse 
x  y  o’o 

component,  therefore  8  -  0  as  shown  in  Fig.  3.  For  uQ  -*•  u^,  (i.e.,  for 
electrons  generated  at  the  tail  of  the  roughness),  8  is  reduced  since  the 
transverse  component  of  electric  field  is  small  there.  Note  that  the  maximum 
value  of  Y88,  located  at  u^/u^  ■  0.4  in  the  case  of  Fig.  3,  is  insensitive  to  $ 
at  fixed  EQ.  This  is  obvious  from 
Eq.  (3): 


d  dpx 

Tt  (Y69)  *  dT 


eE 


m  c 
o 


(6) 


which  approaches  zero  at  distances  sufficiently  far  away  from  the  "bump". 

This  constancy  of  YB9  is  consistent  with  the  notion  that  the  emittance  is 

conserved  during  the  acceleration  process.  The  maximum  value  of  YB8 ,  denoted 

as  Y88  henceforth,  is  independent  of  a  (or  D).  It  is  this  peak  value  YB8 
m  m 

which  will  be  focused  in  the  following,  as  it  provides  a  measure  of  the  beam 

divergence  due  to  the  surface  roughness.  In  fact,  a  simple  argument  would 

show  that  its  use  would  yield  the  upper  bound  for  roughness-induced  emittance. 

The  maximum  value  Y88  depends  on  h,  w  and  E_.  Shown  in  Figs.  4-6  is  the 

m  ° 

dependence  of  Y88m  on  each  of  these  parameters,  with  the  remaining  two 

parameters  being  held  fixed.  The  range  of  these  parameters  is  extended 

arbitrarily  so  as  to  compare  the  numerical  result  with  the  scaling  laws  to  be 

derived  later.  In  Fig.  4,  we  fix  EQ  and  h  on  each  curve.  The  decrease  in 

YB9m  as  w  increases  is  expected  since  enlarging  w  would  reduce  the  transverse 

component  of  the  electric  field.  As  w  <<  h,  the  horizontal  component  of  the 

electric  field  reaches  a  finite  value  near  the  protusion,  so  does  i36  for 

m 

7 


>V-  .s\  \\V* 


small  w.  On  the  curves  in  Fig.  5,  both  E0  and  w  are  held  fixed,  the  maximum 

value  YB0  increases  with  h  for  two  reasons:  (a)  the  transverse  component  of 

m 

the  electric  field  increases  with  h.  (b)  More  importantly,  the  spatial  extent 

of  the  non-uniform  field  increases  with  h,  thus  the  electrons  receive 

transverse  acceleration  for  a  longer  time,  leading  to  a  larger  value  of  YSS^. 

The  increase  of  YSB^  with  EQ,  as  shown  in  Fig.  6,  is  expected  simply  because 

of  the  larger  transverse  acceleration  associated  with  increasing  E  .  Also 

shown  in  the  dashed  curves  in  Figs.  4-6  is  the  value  of  Y38  ,  calculated 

m 

according  to  the  analytic  formula  (15)  given  below.  It  is  seen  from  these 
figures  that  the  analytic  and  numerical  results  have  reasonable  agreement  over 
a  wide  range  of  parameters,  (even  beyond  those  for  normal  operation).  The 
analytic  dependence  of  YB9m  on  h,  w,  EQ  will  become  transparent  in  Section  IV. 
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III.  Spaoe  Charge  Limited  Regime 

It  is  of  interest  to  determine  the  divergence  of  the  electron  beam  in  the 
space  charge  limited  regime,  using  the  same  model  of  the  surface  roughness  as 
in  the  previous  section.  However,  as  remarked  in  the  Introduction,  this  is  a 
rather  difficult  task  if  the  trajectories  of  the  electrons  and  the  resulting 
space  charge  distribution  are  to  be  solved  consistently  over  a  bumpy  surface. 
Analytic  solutions  of  the  Child-Langmuir  type  over  a  protrusion  do  not  seem  to 
have  been  found.  Not  seeking  to  construct  the  Child-Langmuir  solution,  we 
attempt  to  provide  instead  a  rough  estimate  of  the  effects  of  space  charges  on 
the  divergence  of  the  electron  beam,  having  in  mind  the  construction  of  a 
scaling  law  to  compare  with  the  temperature  limited  regime. 

It  is  natural  to  use  conformal  mapping  to  analyze  the  space  charge 
limited  regime:  We  first  use  the  space  charge  limited  electric  field 
distribution  over  a  flat  surface.  We  next  transform  this  electric  field  onto 
the  rough  surface  using  the  conformal  transformation  and  then  calculate  the 
electron  trajectories  subject  to  this  field.  Note  that  several  properties  of 
the  Child-Langmuir  field  are  retained  in  this  procedure:  (a)  The  electric 
field  remains  normal  to  the  conducting  surface,  (b)  The  electric  field  varies 
as  S1^  near  the  surface,  as  in  the  Child-Langmuir  distribution.  Here,  S  is 
the  distance  from  the  rough  surface  in  the  normal  direction,  and  (c)  In  the 
limit  h  +  0,  the  bumpy  surface  approaches  a  flat  surface  and  the  solution 
reduces  to  that  of  a  one-dimensional  Child-Langmuir  flow  whicn  is  truly  a 
consistent  solution. 

Obviously,  the  trajectories  so  constructed  cannot  be  consistent  with  the 
charge  distribution  from  which  they  are  computed.  For  example,  the  electron 
trajectories  in  the  one-dimensional  Child-Langmuir  flow  do  not  intersect, 
whereas  on  a  rough  surface,  the  electron  trajectories  may  do  so.  However,  one 
may  envision  that  the  procedure  described  in  the  preceding  paragraph  as  a 


a 


4 


first  step  in  the  iterative  process:  First,  use  the  one-dimensional  Child- 
Langmuir  electric  field  distribution  over  a  flat  (u,v)  plane  and  transform  it 
onto  the  bumpy  surface  [(x,y)  plane]  as  a  first  approximate  solution. 

Calculate  the  electron  trajectories  over  the  bumpy  surface  subject  to  this 
transformed  Child-Langmuir  field.  [This  is  what  is  done  here.]  The  next 
round  of  iteration  would  be  to  obtain  the  electron  density  from  these 
trajectories,  and  solve  the  Poisson  equation  for  the  electric  field 
distribution  over  the  bumping  surface  due  to  such  a  charge  distribution.  This 
electric  field  distribution  is  a  next  approximation  and  will  be  used  to  re¬ 
compute  the  electron  trajectories  and  the  iteration  process  repeats.  Of 
course,  there  is  no  guarantee  that  such  a  procedure  should  converge  at  all. 
Should  the  iteration  process  converge,  what  we  are  calculating  i3  the  lowest 
order  approximation.  [In  fact,  even  if  this  procedure  can  be  proven 
divergent,  the  result  from  the  lowest  order  iteration  outlined  above  may  still 
yield  useful  information — a  feat  familiar  in  asymptotic  processes.]  To  a3se33 
the  errors  committed  with  the  use  of  this  method,  toward  the  end  of  this 
section,  we  estimate  that  the  bounds  on  emittance  which  we  established 
probably  commit  an  error  on  the  order  30?,  or  less,  over  a  wide  range  of 
parameters. 

With  thi3  caveat,  we  again  consider  a  bumpy  surface  described  by  Eq.  (1  ) 
and  proceed  to  numerically  integrate  the  electron  trajectories  according  to 
(3),  'u),  subject  to  the  initial  conditions  (5).  The  only  modification  is  in 
the  electric  field  (Ex,  Ey)  in  Eq.  (3).  According  to  the  procedure  outlined 
in  the  previous  paragraph,  it  can  be  shown  [Eq.  (B21 )  of  Appendix  B]  that  Ex 


and  E  are  given  by: 


Ex  “  a( -( v  +  b)x  +  u(y  +  b))/[(x<::  +  (y  +  b)2)(u2  +  (v  *  b)*)] 


2  n  1  /2 


Ey  -  a(ux  ♦  (v  ♦  b)(y  ♦  b))/[(x2  ♦  (y  +  b)2)(u2  +  (v  +  b)2)]1/2  (7) 
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under  the  assumption  that  E  -  0  on  the  rough  surface  (space  charge  limited 
regime).  Here,  u  and  v  are  functions  of  x,y  defined  by  the  transformation 
(A1 ),  and 


A  * 


(8) 


In  Eq.  (8), 


Q  »  0.61  * 
*  0.93  * 


r  JD2  1 1 /2 
L(1kA)J 

(p/1  MeV)3A 


(9) 


is  a  dimensionless  parameter  and  the  function  Z(f)  is  the  inverse  of  the 
function  f(Z),  defined  by 

f  -  /  dt  t2(1  +  t*)”1/2.  (10) 

o 

The  current  density  J  in  Eq.  (9)  is  related  to  the  anode-cathode  voltage 

drop  (j)  [See  Appendix  B].  The  approximate  formula  for  Q  given  in  the  last 

expression  of  Eq.  (9)  is  accurate  to  within  5  per  cent  for  all  $  _<  3  MeV. 

1  /2 

[For  p  >  3  MeV,  Q  may  be  approximated  by  (1  ♦  p/.511  MeV)  -  0.8U73,  See 
Appendix  B], 

The  dependence  of  (x  -  xQ)  and  Y89  on  u0  as  the  electrons  reach  the 
anode,  in  the  space  charge  limited  regime,  is  similar  to  Figs.  2,  3-  However, 
since  the  thickness  of  the  space  charge  layer  depends  on  $  (or  J,  or  the 
density)  in  the  space  charge  limited  regime,  the  peak  value  YB8m  depends  on 
four  parameters  h,  w,  <p,  D,  instead  of  three.  Shown  in  Figs.  7-10  are  the 
dependence  of  Y88ffl  on  each  of  these  parameters,  fixing  the  remaining  three  on 
each  curve.  Also  shown  in  these  figures  are  the  analytic  estimates  for  Y80m 
according  to  Eq.  (16)  to  be  derived  in  the  next  section. 
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Note  from  Figs.  7,  8  that  the  dependence  on  h  and  w  are  similar  to  those 
in  Figs.  4,  5.  The  slight  decrease  of  YB0m  with  increasing  $  at  fixed 
Eq  ■  $/D  [Fig.  9]  stems  from  the  fact  that  higher  voltage  means  larger  anode- 
cathode  separation  and,  as  a  result,  the  region  of  low  electric  field  extends 
to  a  larger  distance  from  the  cathode.  Thus,  the  effect  of  the  bump,  near 
which  the  electric  field  is  low,  is  reduced.  This  is  also  reflected  from  the 
curves  on  Fig.  (10,  right).  The  remaining  features  of  these  curves  may 
readily  be  understood  in  a  qualitative  manner  as  in  the  temperature  limited 
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regime  and  will  not  be  repeated  here.  Note  that  YB9m  shown  in  Figs.  7-10  is 
lower  than  the  corresponding  ones  in  the  temperature  limited  regime.  Note 
further  that  the  analytic  formula  yield  reasonable  agreement  with  the 
numerical  results. 

Let  us  now  return  once  more  to  the  error  inherent  in  our  method.  One 
rough  estimate  is  via  the  Poisson  equation 


ifi  +  ifi  ,  SO 


3xe 


3y 


3v‘ 


(n) 


Here,  we  see  that  solving  the  Poisson  equation  in  one  plane  (u,v  plane)  and 
transforming  the  field  in  another  plane  (x,y  plane)  would  commit  an  error  on 
the  order  of  1-J,  where  J  =  | 3(u, v)/3(x, y) |  ,  the  square  of  the  Jacobian.  We 
have  evaluated  J  along  the  trajectory  of  the  electron  which  yields  T39^. 

Fixing  h  *  10  urn,  we  found  that  J  ranges  between  1  and  1.04  for  w  »  JO  h, 
between  1  and  1.16  for  w  -  h  and  between  1  and  0.72  for  w  »  0.1  h.  Thus,  it 
is  plausible  that  the  numerical  value  for  Y38  computed  here,  according  to 
numerical  or  the  analytical  scaling  laws,  are  accurate  to  within  30^  for 
all  h/w  <  10  compared  with  a  consistent  two-dimensional  solution  of  the  Child- 
Langmuir  type. 
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A  ready  assessment  of  the  effects  of  surface  roughness  may  be  obtained  if 


the  analytic  dependence  of  YB8m  is  given  in  terms  of  h,  w,  <j>  and  D.  From  such 
scaling  laws,  a  comparison  between  temperature  and  space  charge  limited  regime 
is  immediate,  at  least  within  the  confines  of  the  present  model.  Such  scaling 
laws  for  both  regimes  are  now  derived. 

First,  Y38  reaches  a  constant  value  at  a  distance  sufficiently  far  away 
from  the  bump,  where  E  -  0  [cf.  Eq.  (6)].  Call  that  distance  L„.  Take 
Ly  »  h.  Since  the  extent  of  the  bump  is  typically  small,  we  write 


Y88  *  88  1  px  1  vx/c. 


The  transverse  velocity  vx  picked  up  by  the  action  of  the  horizontal  electric 
field  over  a  distance  L  is  approximately  given  by 


p  u 

v  »  2  x  Acceleration  x  L  *  E  — s - :  rz  l 

*  *  n  <h2 .  *,2)'/2  x 


where  En  i3  the  normal  component  of  the  electric  field.  In  Eq.  (13),  we  take 
the  directional  cosine  between  the  normal  and  the  x  component  to  be 
h/(h2  +  w2)1/2.  Inserting  (13)  into  (12),  with  Lx  -  h,  we  have 


V33  -  /En  -rr-rirt 

(h  +  w  ) 


where  the  proportionality  constant  in  (14)  is  determined  by  a  "one-point  fit" 
with  the  numerical  solution. 

For  the  temperature  limited  regime,  we  take  E^  »  Eq  -  4>/D  and  fix  the 
proportionality  constant  to  write  (14)  as 

YSem  ■  0.15  /E~  -h-  ■  (rad)  (Temp.  Ltd.)  (15) 

m  0  (h2  ♦  w2)1/4 


W^V.V^V’^WAV.'V 


where  EQ  is  in  MeV/cm,  and  h,w  in  100  urn.  For  the  space  charge  limited 
regime,  we  assign  En  to  be  the  Child-Langmuir  electric  field  evaluated  at  a 
distance  h  from  a  flat  surface,  [cf.  Eqs.  (B13),  ( B 1 6 ) ] .  Inserting  Eq.  (313) 
into  (111)  and  determining  the  proportionality  constant  again  by  a  "one  point 
fit”,  we  have 

1  /2 

Y89m  -  0.079  ■  /4  (jp)2/3(rad)  (Space  Charge  Ltd).  (16) 

(n  +  w  )  ' 

Here  Q  is  given  by  Eq.  (9)  approximately.  The  scaling  laws  (15),  (16)  have 
been  represented  by  the  dotted  curves  in  Figs.  *1-10.  The  agreement  between 
the  analytic  scaling  law  and  the  numerical  calculation  suggests  that  the 
latter  is,  in  fact,  valid  as  long  as  E^  *  S1/3  near  the  rough  surface  for  the 
space  charge  limited  regime.  This  is  the  crucial  assumption  used  in  the 
derivation  of  the  scaling  law  (16).  [All  other  assumptions  employed  in  these 
scaling  laws  have  been  basically  verified  by  noting  that  in  our  analysis  of 
the  temperature  limited  regime,  for  which  the  scaling  law  may  be  taken  as 
reliable,  also  employ  these  assumptions.  The  agreement  between  the  scaling 
law  and  numerical  integration  there  provides  more  confidence  on  their 
validity. ] 

Dividing  Eq.  (15)  by  Eq.  (16),  one  may  arrive  at  the  following  ratio  of 
the  emittances  for  the  temperature  limited  regime  and  the  space  charge  limited 
regime,  for  cathodes  of  tne  same  size  and  type: 

c  (space  charge  limited)  ~  n  r rt -i  1  /6  ,  . 

e  (temp,  limited)  *  5  V 

This  ratio  assumes  a  value  between  0.2  and  0.5  over  a  wide  range  of  parameters 
because  of  the  weak  1/6  power  dependence  in  (17).  This,  perhaps,  is  a  first 
quantitative  assessment  of  the  surface  roughness  effect  in  the  space  charge 
limited  regime.  As  is  evident  here,  it  is  supported  by  extensive  numerical 
calculation. 
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V. 


Discussions: 


Having  derived  an  analytic  scaling  law  for  YB0m  for  an  isolated  roughness 
of  dimension  h,  w,  one  may  readily  deduce  the  statistical  average  <Y&0m>5 

00  00 

<780  >  -  f  f  dhdw  f(h, w)  Y80m  08) 

m  *  J  m 

o  o 

where  f(h,w)dhdw  is  the  probability  of  finding  a  bump  of  height  between  h  and 
h  +  dh,  width  between  w  and  w  +  dw.  Presumably,  f(h,w)  is  to  be  determined 
empirically.  The  value  R^  <Y00m>  would  then  be  expected  to  be  the  upper  bound 
of  emittance  of  an  electron  beam  of  radius  Rb,  due  only  to  the  surface 
roughness.  Simple  estimates  have,  in  fact,  been  obtained  by  assuming 
convenient  forms  of  f(h,w)  and  using  Eqs.  (15)  and  (16)  into  (18). 

The  next  question  is  whether  the  emittance  due  to  surface  roughness  can 
be  measured  in  experiments.  This  depends  on  the  resolving  power  and  on  the 
competing  effects  such  as  nonlinear  space  charge  and  thermal  effects.  Based 
on  the  present  study,  Loschialpo1 ^  and  Shih^  showed  that  surface  roughness  of 
10  um  size  can  be  measured  in  a  two-slit  experiment  using  a  10  kV,  1  A, 
electron  gun  with  slit  width  of  order  1  mil.  The  resolving  power  of  such  a 
slit  limits  the  measurable  brightness  up  to  a  few  times  10^  A/cm2rad2.  Such 
an  experiment  is  being  planned  at  the  Naval  Research  Laboratory.^ 

We  should  stress  once  more  that  our  treatment  of  the  space  charge  limited 
regime  is  only  a  crude  one.  Cur  goal  was  to  provide  a  lowest  order  estimate 
and  to  propose  simple  scaling  laws.  The  usefulness  of  these  formulas  rests  on 
a  direct  comparison  with  an  accurate  solution  (numerical  or  analytical1  of  the 
Chi Id- Langmuir  flow  over  a  bumpy  surface  and,  better  still,  with  experiments. 

A  self-consistent,  two-dimensional  analytic  solution  of  Child-Langmuir  type 
over  a  rough  surface  would  be  highly  desirable.  It  is  perhaps  of  more  than 


academic  interest. 


While  our  calculation  shows  that  operation  in  the  space  charge  limited 


regime  may  reduce  the  deleterious  effect  of  surface  roughness,  the  plasma 
layer  near  the  cathode  surface  may  be  a  source  of  turbulence,  whose  presence 
may  lead  to  a  marked  decrease  of  the  beam  quality. 
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h  -  (b  *  a  ) 


b. 


(A5) 


The  full  width  at  half  maximum,  w,  is 


2b  [ 


3a2  +  2b2  -  2b  (a2  *  b 2)1/2t1/2 
a2  -  2b2  ♦  2b(a2  ♦  b2)1/2 


(A6) 


To  find  a,b  in  terms  of  h,w  from  (A5)  and  (A6)  would  require  numerical 
calculation.  The  following  approximations  obtained  from  (A5)  and  (A6)  may  be 
useful.  For  h/w  >  1 


a  »  h  +  w/3.464 


'■  4/ 


b  *  (w/3.464)[1  +  1/(3.464  h/w)] 


(A7) 


whereas  for  h/w  <<  1 , 


a 


(hw) 


1/2 


(A8) 


b  *  w/2. 


We  have  verified  that  (A7)  is  accurate  to  within  2.5  per  cent  for  all 
h/w  £  1,  whereas  (Ad)  is  accurate  to  within  4  per  cent  if  h/w  <  0.1. 

With  the  branch  appropriately  chosen,  t  is  an  analytic  function  of  z 
which  maps  the  upper  half  5  plane  onto  the  upper  z  plane,  the  latter  being 
bounded  below  by  the  contour  described  in  Eq.  (A2).  Thus  a  uniform 
electrostatic  field  in  the  upper  t  plane  will  be  mapped  into  the  electrostatic 
field  over  the  rough  surface  ^he  temperature  limited  regime.  Such  a 
transformation  is  given,  for  example,  in  Ref.  15.  When  applied  to  the  present 
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model,  the  electric  field  (Ex,  Ey)  in  the  (x,y)  plane  reads 


E  + 
x 


-  iEo  [x  -  i(y  ♦  b)] 

{ [ x  -  i(y  +  b)]2  +  a2}1/2 


(A9) 


which  is  rewritten  as  Eq.  (2)  of  the  main  text.  Here  EQ  is  the  (constant) 
electric  field  far  away  from  the  protrusion.  It  may  be  taken  as  $/D  when  the 
dimensions  h,  w  of  the  bump  are  much  smaller  than  the  anode-cathode  separation 
D. 

From  the  transformation  (A1),  we  may  find  (u,v)  in  terms  of  (x,y).  The 
unit  vector  v  in  the  (u, v)  plane  is  found  to  be  mapped  into  the  vector 


v  -  - 5 - — -r - t — r-rr  |x(-Vx  ♦  uY)  ♦  y(ux  ♦  VY) }  (A10) 

C(u  ♦  V2)(x2  ♦  Y2)]172 

where  V-v+b,Y-y«-b. 
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Appendix  B 


One-Dimensional  Space  Charge  Limited  Flows 


The  one-dimensional  space  charge  limited  flow  (Child-Langmuir  flow) 
serves  as  a  useful  reference  in  the  present  study.  The  relativistic 
formulation  was  given  by  Jory  and  Trivelpiece1 ^  and  is  recapitulated  briefly 
below.  Let  U(s)  be  the  steady  velocity  of  a  one-dimensional  electron  flow  in 
the  s  direction.  The  plane  3  -  0  represents  the  flat  cathode  surface.  Let  J, 
n(s),  and  E(s)  be  the  electron  current  density,  electron  number  density  and 
the  electric  field,  respectively.  The  governing  equations  for  E(s)  and  U(s) 


U  (YU)  -  f-  E(s) 
ds  m 

o 


dE(s)  1_  _J 

ds  e  U( s) ' 

o 


2  2  -1/2 

In  the  force  law  (B1),  Y  «  ( 1  -  U  /c  )  is  the  relativistic  mass  factor  and 

in  writing  the  Poisson  equation  (B2),  we  have  used  the  continuity  equation 
en(s)  U(s)  »  J  =  constant.  Introducing  a  time-like  variable  £,  defined  by 


d£  -  £  -  J  ds/U(  s) 


we  immediately  integrate  (32)  to  obtain 


E  *  —  JE 
Eo 


where  we  assume  E  «  0  at  the  cathode  surface  (s  *  E  -  0). 


We  next  substitute  (35)  into  (31 )  to  obtain 


YU  -  J-  E  -  J  j- 

J  at  me  £ 
o  o 


me 
o  o 


(B6) 


which  leads  to 


cn2Q2/[l 


♦  nVl1/2. 


(B7) 


In  (B7),  both 


n  =  c£/D  (38) 

and 

Q  -  (D/c)  (eJ/2m  ce  )U2  (B9) 

o  o 


are  dimensionless.  Substitution  of  (37)  into  (BU)  yields  3  as  a  function  of 


5: 


s(5)  *  (D/Q)  f  (CQ5/D) 


(BIO) 


where  f(z)  is  defined  as 


f(z) 


dt  t 


1  ,,  .4. 1/2 

0  (1  +  t  ) 


(B11  ) 


whose  inverse. 


which  is  almost  a  constant  for  all  <p  <  3  MeV.  Thus,  to  within  5  pen  cent, 
this  constant  is  0.9.  A  good  approximation  of  Q  then  reads 

Q  *  0.93  (4>/1  MeV)3/4.  (B20) 

Both  Eqs.  (B9)  and  (B20)  are  reproduced  in  Eq.  (9)  of  the  main  text. 

Note  that  all  of  the  well-known  properties  of  the  classical  Child- 
Langmuir  flow  are  recovered  in  the  limit  of  low  voltage.  For  example,  using 
( B 1 6 9  in  (313),  (B14),  we  have  E  ®  s1/3,  $  «  s4^3,  n  «  dE/ds  *  s  2/3  near  the 
cathode  surface  (s  ■*  0). 

In  Section  IV,  we  used  the  one-dimensional  Child-Langmuir  solution  in  the 
C  -  v  +  iv  plane:  [cf.  ( B1 3 ) ] 

2 

„  2m  c 

E  -  v  — QZ(Qv/D)  (B21  ) 

eD 

* 

where  the  unit  vector  v  is  transformed  according  to  Eq.  (A10)  of  Appendix  A. 
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The  peak  value  Y30_  as  a  function  of  the  height  h  in  the  temperature 


limited  regime. 
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The  peak  value  Y86  as  a  function  of  the  electric  field  E_.  in  the 
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temperature  limited  regime. 
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The  peak  value  Y88  as  a  function  of  the  height  h  for  the  space 
in 


charge  limited  regime.  The  analytical  and  numerical  results  are 
practically  indistinguishable  in  this  figure. 


♦  (MeV) 

Fig.  11.  The  dimensionless  functions  Q  and  g(<fi)  defined  in  Eqs.  ( 91 8 )  and 
(B1 9) .  Note  that  g(<i)  is  almost  constant. 


V  V  '^N/  V  ‘.'V/  W\.V_  ■  ...  A  A  A  «*.  /.Aif'v.  A.  A  iVW< 


'iyr  r.'jiyif  »■ 


References 


1.  The  precise  definition  of  electron  beam  quality,  which  is  universally 
adopted  in  the  literature  on  radiation  devices,  is  yet  to  be  established. 
Recently,  Roberson  (Ref.  2)  introduced  one  which  may  be  of  value  to  free 
electron  laser  research.  Earlier  discussions  on  electron  beam  quality, 
beam  brightness,  velocity  spread,  emittance,  etc.,  may  be  found  in  Refs. 
3-5. 

2.  C.  W.  Roberson,  IEEE  J.  Quantum  Electronics  QE-21 ,  860  (1985). 

3.  J.  D.  Lawson,  "The  Physics  of  Charged  Particle  Beams",  (Oxford  University 
Press),  p.  182  (1977). 

4.  C.  Lejeune  and  J.  Aubert,  in  "Applied  Charged  Particle  Optics",  ed.  A. 
Septier  (Academic  Press),  Part  A,  p.  159  (1980). 

5.  J.  Pasour,  R.  Lucey  and  C.  W.  Roberson  in:  Free-Electron  Generators  of 
Coherent  Radiation,  C.  A.  Brau,  S.  F.  Jacobs  and  M.  0.  Seally,  eds., 
[Proc.  Soc.  Photo-Optical  Instrumentation  Eng.  SPIE  Vol.  453,]  p.  328 
(1984). 

6.  S.  E.  Tsimring,  Izv.  Radiofizika,  Vol.  15,  p.  1247  (1972). 


>./ 

V.V 


K 


-« 

:&•: 
f.v.v 

\V.V 

v.  V 
/.W 


V  V  >J 

■xvs 


7.  C.  Kapetanakos,  private  communication  (1985). 

8.  W.  Barletta,  J.K.  Boyd,  A.C.  Paul,  and  D.S.  Prono,  in  Proc.  1984  Free  Electron  Laser  Conf. 
Castelgandolfo,  Italy;  Sept.  1984. 

9.  A.  Shih,  C.  Marrian,  R.  Greene,  private  communication  (1985). 

10.  The  precise  definition  of  electron  beam  emittance  is  given  in  Ref.  4. 

Some  recent  measurements  of  beam  emittance  are  discussed  in  Refs.  9, 

11-13. 

11.  D.  A.  Kirkpatrick,  R.  E.  Shefer  and  G.  Bekefi,  J.  Appl.  Phys.  5[  5011 

(1985). 

12.  A.  C.  Paul,  A.  M.  Sessler,  J.  S.  Wurtele,  G.  J.  Caporaso  and  A.  G.  Cole, 
in  "Free-Electron  Generators  of  Coherent  Radiation",  Eds.  C.  A.  Brau,  S. 

F.  Jacobs  and  M.  0.  Scully,  Bellington,  WA:  SPIE,  p,  108  (1983). 

[SPIE  PROC.  V.  453] 

13-  P.  Loschialpo,  W.  Namkung,  M.  Reiser  and  J.  D,  Lawson,  J.  Appl.  Phys.  57 , 

10  (1985).  Also,  P.  Loschialpo,  private  communication  (1985). 

'4.  Y.  Y.  Lau,  3ull.  Am.  Phys.  Soc.  30,  1=33  (1985). 

15.  F.  3.  Hildebrand,  Advanced  Calculus  for  Applications,  (Prentice  Hall, 

NJ ,  1962),  p.  580. 

16.  H.  R.  Jory  and  A.  W.  Trivelpiece,  J.  Appl.  Phys.  40,  3924  (1969). 


35 


v>V 


m 


y, 


